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Abstract 

Wc construct large families of initial data sets for the vacuum Einstein equations 
with positive cosmological constant which contain exactly Delaunay ends; these are 
non-trivial initial data sets which coincide with those for the Kottlcr Schwarzschild- 
de Sitter metrics in regions of infinite extent. From the purely Ricmannian geometric 
point of view, this produces complete, constant positive scalar curvature metrics 
with exact Delaunay ends which are not globally Delaunay. The ends can be used to 
construct new compact initial data sets via gluing constructions. The construction 
provided applies to more general situations where the asymptotic geometry may 
have non-spherical cross-sections consisting of Einstein metrics with positive scalar 
curvature. 

1 Introduction 

There exists very strong evidence suggesting that we live in a world with strictly 
positive cosmological constant A [38,42]. This leads to a need for a better under- 
standing of the space of solutions of Einstein equations with A > 0. The most gen- 
eral method available for constructing such solutions proceeds by solving a Cauchy 
problem [2, 10, 21]. In view of the general relativistic constraint equations this, sub- 
sequently, requires understanding the corresponding collection of initial data sets. In 
particular one is led to the question of boundary conditions satisfied by the fields. 
When A vanishes a natural set of boundary conditions arises from the obvious model 
solution - the Minkowski space-time. A tempting further restriction is then the re- 
quirement of a well defined and finite total mass, leading to a well understood set of 
asymptotic boundary conditions [1,5,11,34]. When A > the question of asymptotic 
conditions seems to be much less clear cut. One wants to consider a class of space- 
times which includes all solutions of physical interest. Until there is overwhelming 
evidence to the contrary, "physical interest" should carry a notion of "non-singular" . 
The simplest possibility, widely adopted, is to assume that the Cauchy surface 5? is 
a compact manifold without boundary. However, an appealing more general way of 
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ensuring regularity of the initial data is to suppose that (5? ,g) is a complete Rie- 
mannian manifold. One would then like to understand the space of solutions of those 
general relativistic constraint equations with {5? ,g) - complete. 

An interesting class of asymptotic models for such initial data has already been 
explored in the mathematical literature, the time symmetric initial data provided by 
the Delauna$\ metrics [19,20,41]. These describe the family of complete rotationally 
symmetric, conformally flat metrics with constant positive scalar curvature, and are 
in fact well known to general relativists as the time-symmetric slices of the Kottler- 
Schwarzschild-dc Sitter solutions [22,25] (however the connection between these two 
subjects has apparently not been previously noted. In the Riemannian geometric 
context, the Delaunay metrics form the local asymptotic model for isolated singular- 
ities of locally conformally flat constant positive scalar curvature metrics [8,9,24,36] 
(in dimensions n < 5 this also holds in the non-conformally flat setting [26]). The 
known results concerning the existence of complete constant positive scalar curva- 
ture metrics with asymptotically Delaunay ends [7.24.27.30 32.36,37.40] may thus 
be reinterpreted, via their space-time development, as the existence of space-times 
satisfying the Einstein field equations with a positive cosmological constant which 
have asymptotically Kottler-Schwarzschild-dc Sitter ends. 

The object of this work is to point out that every constant positive scalar curva- 
ture (CPSC) asymptotically Delaunay metric is naturally accompanied by a CPSC 
metric with an exactly Delaunay end, and moreover these metrics may be choosen 
to coincide away from the end in question. Such metrics are of interest in general 
relativity for at least four reasons: 

1. They provide, via their maximal development, a large class of space-times sat- 
isfying the Einstein field equations with a positive cosmological constant with 
exactly controlled geometry in the asymptotic regions; in fact the space-time 
development is explicitly known in the domain of dependence of the Delaunay 
regions. 

2. They demonstrate that the special horizon behavior, with alternating cosmo- 
logical and event horizons, which is exhibited by the Kottler-Schwarzschild-de 
Sitter space-time, occurs in large classes of non-stationary solutions. 

3. Any two metrics which carry exactly Delaunay ends with identical mass (De- 
launay) parameters may be glued together using obvious identifications on the 
ends. (A more difficult end-to-end asymptotic gluing theorem of this sort was 
established by Ratzkin [37] , however with exactly Delaunay ends this construc- 
tion is effortless.) Thus Delaunay ends can easily be used as bridges to create 
wormholes, or to make connected sums of initial data sets. Wormholc construc- 
tions are already known to be possible by completely different techniques [16] 
in the setting of a non-positive cosmological constant. Here we provide such a 
construction for positive cosmological constants, with the added bonus of an 
explicit knowledge of the space-time development in the domain of dependence 
of the middle part of the connecting neck, which may be of arbitrary (quantized 
by multiples of the period of the exact Delaunay metric) length. 

4. The asymptotically Delaunay metrics are uniquely characterised by a simple 
geometric criterion [24,26], see Section [2~3l below. 

A natural setting for our considerations is provided by the generalised Kottler 
metrics and generalised Delaunay metrics, as described in Sections l2.2l and l2.4l below. 
Our gluing construction applies in this more general setting. 

1 These are also often called Fowler solutions; see H2.3I for further remarks on the history and choice 
of terminology used here. 
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In an accompanying paper [14], by one of us (PTC) and Erwann Delay, analogous 
constructions are carried out with a negative cosmological constant. With hindsight, 
within the family of Kottler metrics with A G K, the gluing in the current setting 
is the easiest, while that in [14] is the most difficult. This is due to the fact that 
for A > 0, as considered here, one deals with one linearised operator with a one- 
dimensional kernel; in the case A = the kernel is (n + 1) dimensional; while for 
A < one needs to deal with a one-parameter family of operators with (n + 1)- 
dimensional kernels. 

Acknowledgements: DP would like to thank Mihalis Dafermos for first raising 
the question of whether space-times with Kottlcr-Schwarzschild-de Sitter horizon 
behavior exist more generally, and Frank Pacard for a number of illuminating dis- 
cussions. 

2 Kottler— Schwarzschild—de Sitter space and met- 
rics of constant positive scalar curvature with asymp- 
totically Delaunay ends 

In this section we review some results concerning the Kottler-Schwarzschild-de Sitter 
space and CPSC metrics which are asymptotically Delaunay. In order to fix notations 
and conventions we start with some standard facts. 

Recall that initial data for the Einstein field equations with a cosmological con- 
stant A on an n-dimcnsional manifold M consist of a pair (g, K) consisting of a 
Ricmannian metric g on M and a symmetric 2-tcnsor K satisfying the vacuum con- 
straint equations 



where R(g) is the scalar curvature (Ricci scalar) of the metric g. If one considers 
time-symmetric initial data, for which K = 0, then these equations reduce to the 
requirement that g has constant scalar curvature R(g) = 2A. Here we restrict to the 
case where A is positive, and note that the normalization A = "C" -1 ) corresponds to 
R(g) = n(n — 1), the scalar curvature of the standard sphere of radius one in M n+1 . 

2.1 Kottler-Schwarzschild— de Sitter metrics 

The Kottler-Schwarzschild-de Sitter space-time [25] metric in n + 1 dimensions, with 
cosmological constant A > and mass met may be written as 

ds 2 = -Vdt 2 + V-hlr 2 + r 2 h, where V = V(r) = 1 - - ^ , (2.3) 

where £ > is related to the cosmological constant A by the formula 2A = n(n—l)/l 2 , 
while h denotes the standard metric on the unit (n — l)-sphere in R™. To avoid a 
singularity lying at finite distance on the level sets of t we will assume m > 0. 
Equation (|2.3|) provides then a spacetime metric satisfying the Einstein equations 
with cosmological constant A > and with well behaved spacelike hypcrsurfaces 
when one restricts the coordinate r to an interval (rb,r c ) on which V(r) is positive; 
such an interval exists if and only if 



R(g)-(2A+\K\ 2 g -(tr g K) 2 ) = 
Di(K ij - tv g Kg ij ) = 



(2.1) 
(2.2) 




(2.4) 
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When n = 3 this corresponds to the condition that 9m 2 A < 1, and the case of equality 
is referred to as the extreme Kottler-Schwarzschild-de Sitter space-time (for which 
the coordinate expression (|2.3j) is no longer valid). In the limit where A tends to zero 
with to held constant, the space-time metric approaches the Schwarzschild metric 
with mass to, and in the limit where to goes to zero with A held constant the metric 
tends to that of the de Sitter space-time with cosmological constant A. 

The breakdown of the coordinate description above at the horizons r = and 
r = r c can be handled by taking extensions [4,22]: In fact, the Kottler-Schwarzschild- 
de Sitter metric admits an analytic extension (analogous to the Kruskal extension 
of the Schwarzschild metric) as an r-periodic metric on (t,r,6) <G R x M x §" _1 . 
This is most easily seen via the associated conformal Carter-Penrose diagrams [22]. 
The time-symmetric slice t = of the (extended) Kottler-Schwarzschild-de Sitter 
metrics are thus a one-parameter family (parameterized by their mass m) of periodic, 
spherically symmetric, metrics on R x S™ -1 with constant positive scalar curvature 
R= 2A. 

Finally note that, due to the spherical symmetry, each of these metrics is confor- 
mally flat. 



2.1.1 Extreme limit 

It is of some interest to enquire what happens when to — > to, where m denotes the 
values at which equality is achieved in (|2.4j) . In this limit rt, and r c coalesce to a single 
value which we will denote by r. From the space-time point of view the situation is the 
following: recall that the Carter-Penrose diagram for the maximally extended KSdS 
space-times with < to < to is built out of diamond shaped regions corresponding 
to ri, < r < r c , where the Killing vector dt is time- like, and of triangle shaped (either 
upright, or upside-down) regions where dt is spacelike [22]. After passing to the limit 
to — ► to the diamond-shaped regions disappear, and the resulting diagram consists 
of a string of triangles. The Killing vector dt is then spacclike everywhere, except on 
the degenerate horizons r = r^ = r c . 

On the level sets of t a rather different analysis applies, this is discussed in Sec- 
tion [2XU 



2.2 Generalised Kottler metrics 

All the results discussed in Section 12.11 remain valid if to ^ and if the metric h 
in (|2.3j) is an Einstein metric on an (n — 1) -dimensional manifold N 1 '^ 1 with scalar 
curvature equal to (n — l)(n — 2) [6]. We will refer to such metrics as generalised 
Kottler metrics. Note that to = requires (A^ Tl_1 ,/i) to be the unit round metric if 
one does not want r — to be a singularity at finite distance along the level sets of 
t. 



2.3 Delaunay metrics 

The Delaunay metrics, in dimension n > 3, may be defined as the (two parameter) 
family of metrics 

g = u 4 ^ n - 2 \dy 2 + h), (2.5) 

where h is the unit round metric on §" _1 , which are spherically symmetric and have 
constant scalar curvature R(g) = n(n — 1). Thus the functions u = u(y) > must 
satisfy the ODE 

„ in - 2) 2 n(n - 2) ™+2 
u" - ± — —L-u + ^— — '-u^- = 0. (2.6) 
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The two parameters correspond respectively to a minimum value e for u, with 

n — 7 n-2 

<£<£=( )— (2.7) 

n 

(e is called the Delaunay parameter or neck size) and a translation parameter along 
the cylinder. A straightforward ODE analysis (see [32]) shows that all the positive 
solutions are periodic. The degenerate solution with e = corresponds to the round 
metric on a sphere from which two antipodal points have been removed. The solution 
with e = e corresponds to the rescaling of the cylindrical metric so that the scalar 
curvature has the desired value. 

Note that the Delaunay ODE was first studied by Fowler [19,20], however the 
name used here and elsewhere in the literature is inspired from the analogy with the 
Delaunay surfaces: the complete, periodic CMC surfaces of revolution in R 3 [18]. As 
is well known, the analogy between the "conformally flat metrics of constant positive 
scalar curvature" and "complete embedded CMC surfaces of in R 3 " goes far beyond 
this correspondence (see, e.g., [30]). 

Regarding the Delaunay metrics as singular solutions of the Yamabe equation 
on (§™,<7o) one has a number of uniqueness results. Among these are the facts that 
no solution with a single singular point exists, and that any solution with exactly 
two isolated singular points must be conformally equivalent to a Delaunay metric. 
These results can be proved by a generalization of the classical Alexandrov reflection 
argument (the method of moving planes), see [23]. The first general existence result 
for complete conformally flat metrics of constant scalar curvature with asymptotically 
Delaunay ends is due to Schoen [40]. 

Of immediate interest to us is the fact that conformally flat metrics, with constant 
positive scalar curvature, and with an isolated singularity of the conformal factor are 
necessarily asymptotic to a Delaunay metric [24]; in fact, in dimensions n — 3, 4, 5 the 
conformal flatness condition is not needed [26] . Specifically, with respect to spherical 
coordinates about an isolated singularity of the conformal factor, there is a half- 
Dclaunay metric which g converges to, exponentially fast in r, along with all of 
its derivatives. This fact is used in [27,31,32,36,37] where complete, constant scalar 
curvature metrics, conformal to the round metric on S\{pi, . . . ,Pk} were studied and 
constructed. (This is one instance of the more general "singular Yamabe problem" .) 

By the uniqueness of solutions to ODEs, or otherwise, we have: 

Proposition 2.1 The time symmetric initial data sets for Kottler-Schwarzschild- 
de Sitter space in spatial dimension n with A = n ( Tt ~ 1 ) ) are precisely the Delaunay 
metrics with constant positive scalar curvature R = n(n — 1). 

This correspondence continues to hold for any choice of positive cosmological 
constant A provided that one homothetically rescales the Delaunay metrics so that 
R = 2A. 

Comparing (|2.3p and (|2.5|) we find 

r = u^- , r^- = V~ 1/2 , (2.8) 
dr 

which allows us to determine y as a function of r on any interval of r's on which V 
has no zeros. 

2.3.1 Extreme limit 

Let m and f be as in Section 12.1.11 and suppose that < m < m, denote by r* £ 
(rb, r c ) the value at which the maximum value V* of V is attained, shifting y by a 
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constant we can assume that the corresponding value = y(r#) of the y coordinate 
in (|2.6|) is zero. We have r, -> r and K. — ► as m — ► to, and it clearly follows from 
(|2.8j) that the correspondence y <-> r breaks down in the limit. This singular behavior 
with respect to the r coordinate is of course resolved by the coordinate y of (|2.5p . 
A somewhat more explicit way of seeing this is to replace r by a new coordinate w 
through the formula 

r = r* + vKw 



which scales up the interval r G (r&, r c ) to w G ((r& — ?**)/v^4, (f c — 
Equations (|2.8p become 

r* + v^' = , t L = 7 ^ : r- , (2.9) 



which are regular in the limit m — > rh. In the new coordinates we have 

-^pdr 2 + r 2 h = ^dw 2 + (r* + \/V^w) 2 h > dw 2 + r 2 h as m — > m , 

with the limit being uniform over compact sets of the iu coordinate. This shows in 
which sense the space sections of the KSdS metrics approach a cylindrical geometry 
in the extreme limit. It should, however, be borne in mind that the space-time picture 
of Section 12.1.11 is rather different. 



2.4 Generalised Delaunay metrics 

Similarly to Section l2~2l the analysis presented at the beginning of Section l2.3l rcmains 
valid when the parameter e is positive and if the metric h in (|2.5p is an Einstein metric 
on an (n—1) -dimensional manifold N"^ 1 with scalar curvature equal to (n— l)(n— 2). 
We refer to the resulting metrics as generalised Delaunay metrics. As before, e = 
requires (N n ~ 1 1 h ) to be the unit round sphere if one wants to avoid a singularity at 
the set u(y) = 0. 



2.5 Complete metrics with constant positive scalar curvature 
and asymptotically Delaunay ends 

Conformal gluing constructions for constant scalar curvature metrics g = u 4 ^ n - 2 ^g 
with R(g) = R(g) = n(n — 1) have given rise to a wide variety of such metrics 
with asymptotically Delaunay ends. The linearisation of this equation about a so- 
lution leads to the operator L g = A g + n. The key nondegeneracy assumption of 
any conformal gluing construction is that L g is surjective when acting on appro- 
priately defined function spaces. The Delaunay metrics themselves are nondegen- 
cratc in this sense [32], moreover the solutions constructed by Mazzeo-Pacard on 
§ \ {pi, . . . ,Pk} [27] are non-degenerate. On the other hand, the standard metric on 
the n-sphere, (§, go), is degenerate due to the fact that the restrictions of the lin- 
ear functions in R™ +1 span an (n + l)-dimcnsional co- kernel of L go . In addition to 
the original construction of Schocn [40], the constructions of [31] and [37] use non- 
degenerate solutions as building blocks to produce new non-degenerate solutions. 

All of the constructions alluded to above are in the setting where the metrics are 
locally conformally flat everywhere. This is clearly not necessary. A general conformal 
gluing theorem was established by Byde [7]: 

Theorem 2.2 (Byde [7]) Let {M,g) be a compact Riemannian manifold, possibly 
with boundary, of constant scalar curvature n(n — 1), which is non- degenerate in the 
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sense described above, and let xq € int(M) be a point in a neighborhood of which 
g is conformally flat. Then there is a constant po and a one parameter family of 
complete metrics g p on M\{xq} defined for p G (0, po), conformal to g, with constant 
scalar curvature n(n— 1). Moreover, each g p is asymptotically Delaunay and g p — > g 
uniformly on compact sets in M\ {xq} as p — > 0. 

This result is exactly analogous to results for constant mean curvature surfaces 
established in [28,29]. Byde goes further and shows how one can also glue asymptot- 
ically Delaunay ends onto non-compact, non-degenerate solutions (though without 
the uniform convergence to the original metric away from the gluing locus). 

Note that, in light of Proposition 12.11 all of these results, and others, on the 
existence of CPSC metrics with asymptotically Delaunay ends have an immediate 
reinterpretation, after considering the maximal development of the initial data set, 
as statements regarding the existence of space-times with asymptotically Kottler- 
Schwarzschild-de Sitter ends. 

3 Perturbation to exactly Delaunay, or generalised 
Delaunay ends 

The gluing construction of Corvino-Schoen [17,39] (compare [13]) generalises to the 
positive cosmological constant setting as follows: 

Theorem 3.1 Let iV™^ 1 be compact, let (M, g) satisfy R(g) = n(n— 1), and suppose 
that M contains an end E m [0, oo) x TV" -1 on which g is asymptotic to a generalised 
Delaunay metric g = g £ , with < e < ( 2 ^)^^; together with derivatives up to 
order four. Then for every S > there is an e' satisfying \e — e'\ < 5 and a metric 
g' with R(g') = n{n — 1), which differs from g only far away on E, and which is a 
generalised Delaunay metric with Delaunay parameter e' on the complement E' of a 
compact subset of E. 

Remark 3.2 By taking the maximal, globally hyperbolic, space-time development 
of the time-symmetric initial data set (M, </), we obtain a solution of the vacuum 
Einstein equations with cosmological constant A = n(n — l)/2 such that the metric 
on the domain of dependence of the end E 1 is isometric to a subset of the Kottler- 
Schwarzschild-dc Sitter space-time. 

Proof: Let g asymptote to a generalised Delaunay metric g e on R x N n ^ 1 . We can 
write g e as 

<? e = dx 2 + e 2f &h , (3.1) 

where h is an Einstein metric on iV n_1 , normalised as described above. 

Consider a connected component of the set on which V > 0, where V is the 
function appearing in (|2.3p for the metric g E . It follows from (|2.3[) that \/V is the 
normal component of the Killing vector d t on the level sets of t. From the general 
results in [33] it follows that any such function, for a static space-time, solves equation 
(13. 2|) below. It further follows from the analysis in [12] that VV can be smoothly 
continued to a real-analytic function on M, which we call N, by changing signs across 
the zero level sets of V. Furthermore, both N and V = ±N 2 are functions of x only 
in the representation (|3.ip of g. 

Let T = T(e) be the period off, and let ^ = [iT + a, iT + T + a] x TV™- 1 , where 
a will be chosen below. 
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Let (R x N n 1 ,g e ') be a generalised Delaunay metric with parameter e' near e. 
\e-e'\ < S. 

Let g £ i be a metric on f2j obtained by interpolating between g and <? e < using any 
i independent cut-off function smoothly varying from zero to one. The cut-off should 
be supported away from the end-points of the interval [i, i + T]. 

To achieve constant scalar curvature we will, first, correct the metric g e > to a 
new metric g E > using the operator L of [13], as restricted to time-symmetric data, so 
that Y = there. The correction will be of the order of the perturbation introduced, 
namely 0(\e — e'\). This will, however, not quite solve the problem because the 
operator L at g = g e has a cokernel, which consists of functions solving the "static 
KIDs equation" : 

DiDjN = NR l0 + A g N gij . (3.2) 
We thus have to understand the space of solutions of (|3.2|) : 

Lemma 3.3 Let g = g be a generalised Delaunay metric as in \3.1\) . 

1. If (N n ,h) is the round sphere and if m = 0, then the space of solutions of 
\3. 2]) is (n+ 1) -dimensional, spanned by N together with functions of the form 
e$ a 1 , where f is as in LS. 1\) and a 1 is the restriction of the Euclidean coordinate 
x l to § n_1 under the standard embedding S n_1 R n . 

2. Otherwise all solutions of \3. 2\) are proportional to N . 

Proof: Let v A be coordinates on the level sets of a;; we have T X AB = —f'hAB, 
Y% 3 = 0, T^Q = f'Sg, and T BC = tg C , where the f g C 's are the Christoffel symbols 
of h. Since N depends only upon x, and satisfies (|3.2j) . we immediately find 

R xA = 

away from the zero-set of V; by continuity this holds everywhere (this conclusion 
could also have been reached directly from the warped product structure of g) . But 
then (|3.2[) gives 

= D X D A N = e f d x (e- f d A N) , 

hence 

N(x, v A ) = e 1[x) P{v A ) + M(x) , (3.3) 

for some functions P = P(v A ), M = M(x). 
Set 

a := N/N , 

then a is smooth away from the zero-level sets of N. Since both N and N satisfy 
(|3.2|) one finds that a is a solution of the equation 

ND i D j a + D i aD j N + D j aD i N = 0. (3.4) 

From DaN = we obtain 

D A D B a = 0. (3.5) 

Let Xab = — F\b — f^AB be the second fundamental form of the level sets of x, 
(|3.5p can be rewritten as 

§>A®bol + a x X AB = ® A 2>Ba + Oxf^hAB = , (3.6) 
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where is the covariant derivative operator of the metric h, and a x = d x a. Applying 
$) B to (|3.6|) and commuting derivatives one obtains (recall that the Ricci tensor of 
h equals (n — 2)h) 



Contracting A with B in (|3.6[) we find ip = 3>b 
conclude that there exists a constant C such that 



J B @ B a + (n-2)a- tp) = 0. 

B a/(n — 1), which allows us to 

LStant C such that 
B (a + C) = -(n- l)(a + C) . (3.7) 



Suppose, first, that (A n_1 ,/i) is the unit round sphere. Equation (|3.7[) shows that 
solutions of (|3.6[) are linear combinations of the constant function a = 1 and of the 
functions a 1 = x i \§n-i, where x l is a canonical coordinate in R n , with S™" 1 being 
embedded in R™ in the obvious way. Hence, there exist functions A^(ai), fJ, = 0, . . . , n, 
such that 

a(a;,w^) = A M (x)a M (u^) . (3.8) 

But then (|3.3p implies \i(x) = e^N^ 1 '^ for some constants Ai, without however 
imposing any constraints on Xo(x) which remains undetermined so far. 

On the other- hand, if (iV n-1 , h) is not the unit round sphere, then by a theorem 
of Obata [35] the function a does not depend upon v A , so that (|3.8|) again holds with 
a 1 = and a = 1. 

Inserting (|3.8[) into (|3.4p with ij = one finds 

^{iV 2 ^)^ [Ao(x) + N- 1 (x)e f ^X l a t (v A )] } = . (3.9) 

with p.ip to conclude 

(3.10) 



In order to analyse this equation, it is useful to compare 
that 



dr _ 1 



1 



We then have TV 2 ^ = ±A> 3 <9 r and 



±iv 2 a r 



Ao + A^^^/^a 4 ^) 



N A d r 



2m 



A + 



2m 



=A ia > A ) 



iv 3 9 r A + 1 - ^Ua > A ) • 



So p.9p will hold if and only if this is a function which depends at most upon v A . 
Hence Ao is a constant and A^ = unless m = 0, in which case the Ai's are arbitrary, 
as desired. □ 



Returning to the proof of Theorem 13. 1[ the metric g e i is obtained by solving the 
equation 

(ii(SfeO - n(n - 1)) | n _ € (Iml) 1 

using the implicit function theorem, compare [13, Theorem 5.9]; this can be done 
on Qi for all i large enough. As already mentioned, the perturbation introduced is 
0(\e — e'\). In view of Lemma I3~3l the obstruction to solving the problem is thus the 
vanishing of 

N(R(g e ,)-R(g e ,))dM, (3.11) 
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where N is the static KID associated with the generalised Delaunay metric g e . We 
need the following identity, from [15]: 

Vdct^ N(R g - R b ) = d, (u'(iV)) + Vdrtg (p + Q), (3.12) 

where 

D*(JV) := 2Vdit^ (Ng i V t gN@ j g H + D^Ng^ k e jk ) , (3.13) 

p := (-NRic(b)ij + &&jN - AbNbi^gtogf'ekt , (3-14) 
Q := N(gv - 6« + 3 ifc 5 ^e w )Ric(6) y + Q' . (3.15) 

Brackets over a symbol denote anti-symmetrisation, with an appropriate numerical 
factor (1/2 in the case of two indices). Here Q 1 denotes an expression which is bilinear 
in 

e = eijdx l dx° :— (gij — bij)dx t dx J , 

and in ^e^, where @ denotes now the covariant derivative operator of the metric b, 
linear in N, dN and HessTV, with coefficients which are constants in an ON frame for 
b. The idea behind this calculation is to collect all terms in R g that contain second 
derivatives of the metric in diW; in what remains one collects in p the terms which 
are linear in e^, while the remaining terms are collected in Q; one should note that 
the first term at the right-hand-side of (|3.15[) docs indeed not contain any terms 
linear in when Taylor expanded at gij = bij. Note that p vanishes when b = g e 
by choice of N. So the integrand is quadratic in ey, up to terms 0(|e — s'\)eij, and 
up to the divergence which produces a boundary term 

WdSi . 

For our next lemma it is convenient to write two generalised Delaunay metrics g 
and b as 

dr 2 9 ° dr 2 « . 

9=x2+r 2 h, b=^+r 2 h. (3.16) 

We claim: 

Lemma 3.4 Let g and b be two generalised Delaunay metrics with mass parameters 
m and itiq. Let r be such that N(r) ^ and N(r) 7^ 0. If {r} x N"^ 1 is positively 
oriented, then 

[ U i dSi = 2u n - 1 (n-l)NN- 1 (m-m Q ) . (3.17) 

•/{rjxAT™- 1 

where w n ~i is the volume of N^ 1 . 

Proof: Let us denote by r*- ft the Christoffel symbols of the metric b. We have 
T r rr = -d r N/N, Y r AB = ~rN 2 h AB , F^ c = f^ c (where, as before, the F^ c 's are 
the Christoffel symbols of the metric h), Tf B = r~ 1 5 B l , while the remaining F's 
vanish. It holds that e = (N~ 2 - N~ 2 )dr 2 , fr om which one easily finds 



U r = 2(n - l)NN~ 1 (m~m )Vdeth\ (3.18) 
and the result follows by integration. □ 
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We are ready to show that one can choose e' — equivalently m! — so that the 
obstruction vanishes. So we consider the integral (|3.11|) . We wish to use (|3.12[) with 
g = g £ i . Note that the integration in (|3.11|) is taken with respect to the measure dfi g , 
while (|3.12p involves d\i g . The difference between the two volume integrals comes 
thus with a prcfactor 0{\e — e'\), and produces an error term which is 0((e — e') 2 ) : 

J N[R{g s .)-R{g e ,))d^ 

= j N(R&) - R(g e >j)d(i g + 0((e - e') 2 ) 

= [ U l dSi- [ WdSi + 0((e-e') 2 ) . (3.19) 

We now choose a in the definition of fli so that the number 

A:=4w n _ 1 (n-l)iV| ;c=(r 

does not vanish; note that A equals, up to 0(\e — s'\), the number in front of (to — too) 
in (|3.17[) . Since g approaches g together with its first derivatives as i goes to infinity, 
the first integral in the last line of (|3.19j) is o(l), where o(l) tends to zero as i 
tends to infinity. By Lemma T3. 41 the second integral in the last line of (|3.19|) equals 
A(m' — m) = 0{\e — e'|), where to' is the mass parameter of 7 e / while m is that of g. 
We infer that 

J N(R{g e .) - R(g s '))dfi g = A(m' - m) + o(l) + 0((e - e'f) . 

Clearly this can be made positive or negative when i is large enough by choosing e' 
appropriately; by continuity there exists an s' which makes the integral vanish, and 
the result is proved. □ 



4 Concluding remarks 

Our work leads naturally to the following questions: 

1. In light of the results in [26] one should expect that, at least for n = 3, 4, 5, a con- 
struction similar to Byde's [7] could be carried out without any assumption of 
conformal flatness in a neighborhood of the omitted point (which forms the end 
of the resulting complete metric). Moreover, with or without the conformally 
flat condition, it should be straightforward to iterate Byde's construction to 
produce any number of asymptotically Dclaunay ends. If one could add asymp- 
totically Dclaunay ends at any chosen set of points in any positive constant 
scalar curvature manifold, then our analysis here could then be used to replace 
them with exactly Delaunay ends. (Furthermore, all of this should be doable 
via a local deformation of the metric near points without static KIDs, using the 
techniques of [3, 13, 16].) Alternatively, can one generically deform a constant 
positive scalar curvature metric, keeping the scalar curvature fixed, to a metric 
which is conformally flat near a set of prescribed points? Proposition 4.1 of [41] 
could perhaps be used as an intermediate step here. 

2. In Byde's construction, or in a variation thereof as just suggested, can one 
ensure that the range of masses of the resulting Delaunay ends covers an interval 
of the form (0, e), for some e > 0? This is a natural condition which has appeared 
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elsewhere as a necessary hypothesis (see e.g. [37]). Now, it is clear that the 
masses of our exactly Delaunay ends are continuous functions of the initial 
mass. Given any two points with the associated families of exactly Delaunay 
ends, one could then always adjust the masses to be the same, ensuring that 
the ends can be glued together. 

3. We did not carry out the gluing in situations when the metric g approaches a 
cylindrical metric dy 2 + h along the asymptotic end; such metrics arise in black- 
hole space-times with degenerate horizons. This deserves further attention. 

4. It would be of interest to extend the current gluings to general rclativistic initial 
data with non-vanishing extrinsic curvturc. 
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